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Abstract Instances of local damage in timber such as knots, 
decay, and cracks can be translated into a reduction of 
service life due to mechanical and environmental loadings. 
In wood construction, it is very important to evaluate the 
weakest location and to detect damage at the earliest pos¬ 
sible stage to avoid future catastrophic failure. In this study, 
modal testing was used on wood beams to generate the first 
two mode shapes. A novel statistical algorithm was pro¬ 
posed to extract a damage indicator by computing mode 
shapes of vibration testing before and after damage in tim¬ 
bers. The different damage severities, damage locations, 
and damage counts were simulated by removing mass from 
intact beams to verify the algorithm. The results showed 
that the proposed statistical algorithm is effective and suit¬ 
able for the designed damage scenarios. It is reliable for the 
detection and location of local damage of different severi¬ 
ties, location, and number. The peak values of the damage 
indicators computed from the first two mode shapes were 
sensitive to different damage severities and locations. They 
were also reliable for the detection of multiple cases of 
damage. 

Key words Damage detection • Mode shape • Timber • 
Nondestructive test 


Introduction 

The average modulus of elasticity (MOE) of wood, which is 
determined by global nondestructive test techniques, is 
commonly used to design its usage by considering design 
safety factors. The MOEs at various locations are mostly 
identical for a uniform wood beam without any defects. 
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However, most solid wood beams contain natural defects 
such as knots so that its MOEs at various locations vary. 1 
Wood beams are also susceptible to fungal decay and crack 
development during their service life due to mechanical and 
environmental loads. Such forms of local damage are 
mostly described as failure causes of reduced service life. If 
the weakest portion can be determined using nondestruc¬ 
tive tests, the safety design factors can be optimized to 
improve the usage of wood products. It is also very impor¬ 
tant to monitor the weakest location and detect damage at 
the earliest possible stage to avoid the future catastrophic 
fracture in most engineering applications. 

Over the past 30 years, damage detection in structures 
has received considerable attention from the civil, aero¬ 
space, mechanical engineering, and wood engineering com¬ 
munities. Commonly used damage detection methods are 
either visual or localized experimental methods, such as 
acoustic, ultrasonic, magnetic field methods, radiograph, 
eddy current, and thermal field methods. Vibration-based 
methods in wood engineering were studied by some 
researchers in the past decade to detect dynamic proper- 
ties “ and local damage. ’ “ The basis of vibration-based 
damage detection methods is that any changes in the physi¬ 
cal properties of the wood beam, i.e., stiffness, mass and/or 
damping, will in turn alter its dynamic characteristics (reso¬ 
nant frequencies, modal damping, and mode shapes). Yang 
et al. reported the possibility of using the shape of the 
flexural vibration and its curvature of the first mode to 
detect the damage. However, it is still necessary to study the 
damage pattern because it is vital for a damage detection 
algorithm that the location and severity of defects and/or 
damage should be automatically identified by certain dam¬ 
age indicators. In addition, a damage-detecting algorithm 
must be suitable for multiple defects and/or forms of 
damage. 

In this study, a novel statistical algorithm is proposed to 
compare mode shapes of vibration testing for damage de¬ 
tection and location in free-free supported wood beams. 
The proposed method successfully verified the different 
simulated scenarios of different damage locations, different 
damage severities, and multiple damage points. 
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Table 1 . Simulated damage scenarios 


Scenarios 

Case 1 

Case 2 

Case 3 

Case 4 

Case 5 

Case 6 

Case 7 

Case 8 

Case 9 

Case 9 

w (mm) 

3 

3 

3 

3 

3 

3 

3 

3 

3 

3 

1/L 

40% 

40% 

40% 

40% 

25% 

40% 

55% 

70% 

40% 

70% 

Impact point 

16-17 

16-17 

16-17 

16-17 

7-8 

16-17 

25-26 

34-35 

16-17 

34-35 

h/H 

5% 

10% 

20% 

40% 

20% 

20% 

20% 

20% 

20% 

20% 


Beam 1 for cases 1-4; beams 2-5 for cases 5-8; beam 6 for case 9 


FFT analyzer 


Laptop 


Amplifier 


Accelerometer 

Beam 



Fig. 1 . Diagram of the detecting system 


Materials and methods 

Six uniform and clear beams without any defects after being 
visually graded were cut from a log of white oak ( Quercus 
rubra) and their moisture contents were conditioned to 
around 12%. The beam dimensions were 1200 x 25.4 x 
50.8mm (length x width x thickness). The damage scenarios 
were simulated by removing mass from the intact beams 
and the different damage scenarios are listed in Table 1. 

The schematic diagram of the experimental appartus is 
shown in Fig. 1. The measuring system consisted of an im¬ 
pact hammer, an accelerometer, two preamplifiers, and a 
dual-channel fast Fourier transform (FFT) analyzer. Each 
beam was tested before and after damage. The beam was 
free-free supported, 14 and the part between the supporting 
points was subdivided into 45 small equal block elements by 
44 center-symmetric impacting points along the algebraic 
center point of the beam length. The beam was impacted by 
the hammer, point by point, and the resonant signal was 
picked up by an accelerometer, which was fixed at one of 
the free ends. Both the output signal from the hammer and 
the accelerometer were amplified by preamplifiers and 
recorded simultaneously using the dual-channel FFT ana¬ 
lyzer. Then the data were dumped from FFT analyzer to 
personal computer for postprocessing. The postprocessing 
software was developed in MatLab program language for 
damage identification. 


Damage-locating algorithms 

The natural frequencies and mode shapes were obtained 
from the recorded time domain data by computing the 
frequency response function (FRF) for each measuring 
point. The Fourier transforms of excitation (input) and re¬ 
sponse (output) signals can be expressed by the following 
equations: 


U(co) = J u(t)e- ,m dt 



V(co) = J v(t)e-“°'dt 



where u(t) is the excitation signal outputting from the sen¬ 
sor hidden inside the hammer; v(t ) is the response signal 
recorded by the accelerometer; U(co) is the Fourier trans¬ 
form of u(t ); and V(co) is the Fourier transform of v(t). 

The frequency function [FRF or H(co)] of a system was 
computed as follows: 



Equation 3 can be rewritten as the following equation: 


H(co) = H 1 (co) 


G uv (cq) 
G uu (co ) 



where G uv (<n) = U*(co)V((0) is the cross-spectrum between 
u(t) and v(t)\ and G vv (co ) = U*(cq)U(cd) is the autospectrum 
of u(t). 

The amplitude of the peak (/y) in the FRF at frequency 
cOjj for measuring point j is the modal coordinate for measur¬ 
ing point j at mode shape i. Here the first two mode shapes 
were studied so that i- 1,2 and j - 1,2,..., 44. The mode 
shape coordinates for all measuring points were connected 
to form the mode shape i. The difference of the mode 
shapes after and before damage was expressed as follows: 


— 77 aam _ 77 
r ij ~ r ij F ij 


’dam 


int 



where F^ f is defined as the damage indicator, which is the 
difference of the mode shapes i at the y-th measuring point 
before and after damage; Ff* m is the mode shape coordi¬ 
nates after damage at the y-th measuring point for mode 
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shape i; F is the mode shape coordinates before damage 
at the y-th measuring point for mode shape i. Discrete 
Laplacian transform expressed by Eq. 6 was then operated 
on the difference after the above-mentioned subtraction. 



Considering the difficulty in discerning whether the damage 
indicator £ tj indicates damage at the y-th location by compar¬ 
ing the damage indicator for one element with another, the 
damage indicator values (£ i} ) were assumed to be random 
variables and were evaluated using a standard normal statis¬ 
tical model. Standard normal damage indicator values are 
denoted as follows: 


where DL, is the standard normal indicator value of the 
mode shape i at the y-th location; and /r, and o; are mean and 
standard deviation of for all measuring points at the 
mode shape i, respectively. The element was classified as 
damaged if its standard normal indicator value was greater 
than a threshold, which was selected judgementally. 


Results and discussion 

The first two natural frequencies of six beams before and 
after damage are listed in Table 2. It can be clearly observed 
that the natural frequencies after the removal of some mass 
are smaller than those of intact beams except for the first 
frequencies of cases 1 and 2. This is because the stiffness of 
the beam is reduced after removing mass. Some researchers 
reported that the natural frequency is a potential nonde¬ 
structive indicator to evaluate and grade timbers. 14 As 
shown in Table 2, the first natural frequencies of cases 1 and 
2 showed no difference when compared with intact beam. 
These results indicate that the first natural frequencies are 
not suitable for detecting low-level damage or disturbance, 
while to date there is no known means of using the natural 
frequencies to locate local damage. 


Table 2. The first two natural frequencies before and after damage 

Scenarios 

First frequency 

Second frequency 

Before 

damage 

After 

damage 

Before 

damage 

After 

damage 

Case 1 


113.0 


301.9 

Case 2 


112.8 


299.1 

Case 3 

113.0 

109.5 

306.1 

287.1 

Case 4 


104.4 


280.3 

Case 5 

119.1 

115.1 

323.9 

314.9 

Case 6 

117.1 

111.0 

319.2 

294.5 

Case 7 

146.6 

135.7 

363.5 

357.9 

Case 8 

124.1 

118.0 

334.3 

316.0 

Case 9 

113.3 

109.2 

305.9 

297.8 


Frequencies given in hertz 


The different damage severities (cases 1-4) were simu¬ 
lated to study the effects of the damage severity on the 
sensitivity of the proposed damage indicator. Damage indi¬ 
cators computed from the first two mode shapes versus 
impacting points under different damage severities {h!H\ 
see Fig. 1) are shown in Figs. 2 and 3, respectively. It can be 
observed that the damage can be located by the damage 
indicators computed from the first mode shape or the 
second mode shape. It is reasonable to conclude that 
the proposed algorithm is sensitive to damage severity. The 
damage indicator increases when the damage severity in¬ 
creases. The peak values of the damage indicator computed 
from the first mode shape under different damage severities 
(h/H: 5%, 10%, 20%, and 40%) were 2.04, 2.2, 2.29, and 
2.63. The peak values of the damage indicator computed 
from the second mode shape were 2.26, 2.43, 2.44, and 2.84. 
Although the above-mentioned information should not be 
construed to assume that the damage severity can be evalu- 



Fig. 2. Damage indicators (first mode) versus impacting points under 
different damage severities of 5%-40% 



Fig. 3. Damage indicators (second mode) versus impacting points 
under different damage severities of 5%-40% 
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Fig. 4. Damage indicators (first mode) versus impacting points under 
different damage locations (x/L) of 25%-70% 



Impact point 

Fig. 5. Damage indicators (second mode) versus impacting points 
under different damage locations (x/L) of 25%-70% 


ated quantitatively and exactly, they can be used to identify 
the existence of damage and its location. Further studies 
should be conducted to improve the algorithm so that the 
damage severity and the remaining useful life can be esti¬ 
mated in structures. 

It is well known that the boundary conditions and the 
mounting of sensors should be carefully considered when 
designing a vibration test. In this study, different damage 
locations (cases 5-8) were simulated to study the effects of 
different damage locations on the effectiveness of the pro¬ 
posed damage indicator. Damage indicators of the first two 
modes versus impacting points under different damage lo¬ 
cations (x/L; see Fig. 1) are shown in Figs. 4 and 5, respec¬ 
tively. Generally, it can be reasonably concluded from Figs. 
4 and 5 that the simulated damage at different locations 
along the beam were correctly detected and located by the 
proposed indicator (DI, ;/ ). The peaks of damage indicators 



Fig. 6. Damage indicators versus impacting points with multiple dam¬ 
age points 

coincided well with the designed damage locations, and the 
existence of damage and its location can be readily distin¬ 
guished by DI /; . The damage indicators tended to become 
bigger under the first mode, while the damage indicators 
became smaller under the second mode when x/L was 40% 
and 55%. It was thought that this interesting phenomenon 
was due to the effect of the artificial damage on the local 
stiffness of the beam and the nature of the mode shapes. 
The artificial damage in this study decreased the natural 
frequencies of the beam as shown in Table 2 due to the 
reduction of the local stiffness of the beam. Meanwhile, the 
theoretical mode shapes follow a certain sinusoidal function 
for a free-free supported beam. 15 The maximum deforma¬ 
tion of the first mode shape approaches 50% (x/L) and the 
deformation tends to be zero at both ends in a sinusoidal 
pattern. DI is bigger when xlL is 40% or 55% than that 
when xlL is 25% or 70% as shown in Fig. 4. The reason 
is that the deformation is bigger when xlL is 40% or 55%. 
The maximum deformation of the second mode shape ap¬ 
proaches 25% and 75% (x/L) and its deformation is zero at 
50% (x/L) and at both ends in a sinusoidal pattern. A simi¬ 
lar interpretation can be applied to explain why DI is bigger 
when x/L is 25% or 75% in Fig. 5. The damage indicator 
extracted from the mode shape is meaningful only when the 
deformation and/or deflection are not zero. Therefore, at 
least the first two mode shapes are suggested to be used 
together for damage detection to avoid the effect of the 
mode shape. However, it was found that it was difficult to 
generate the higher mode shapes such as the third mode 
shapes because of the excitation source or noise in the 
response signal. 

Two damage counts (case 9) were also simulated to 
verify the effectiveness of the proposed damage indicator. 
Damage indicators computed from the first two mode 
shapes versus impacting points under two damage counts 
were correctly detected as shown in Fig. 6. It can be ob¬ 
served that only the damage indicators of several impact 
points close to the damage are relevant to the disturbance. 
It is readily understood that there is no interaction between 
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the multiple damage points with certain distances because 
the mode shape-based damage indicator is a local indicator 
and not a global indicator. Therefore, the mode shape- 
based algorithm is suitable for detection of cases of multiple 
damage areas. 


Conclusions 

A statistical algorithm was proposed for damage detection 
in timbers. The different severities of damage, damage 
locations, and damage counts were simulated by removing 
mass from intact beams to verify the algorithm. The results 
showed that the proposed algorithm is effective and suitable 
for the designed damage scenarios. The algorithm is reliable 
for the detection of local damage under different severities, 
locations, and counts. 
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